Measurements of the Dgj resonance properties 
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We report measurements of the properties of the _D^j (2317) and -D^j(2457) resonances produced in 
continuum e^e~ annihilation near y/s = 10.6 GeV. The analysis is based on an 86.9 fb~^ data sample 
collected with the Belle detector at KEKB. We determine the masses to be M (Z)+^(2317)) = 2317.2± 
0.5(stat)±0.9(syst)MeV/c2 and M(D+_;(2457)) = 2456.5± 1.3(stat)± 1.3(syst) MeV/c^ We observe 
the radiative decay mode D^j(2457) — > -D^7 and the dipion decay mode _D^j(2457) — » D^7r"'"7r~, 
and determine their branching fractions. No corresponding decays are observed for the Dsj (2317) 
state. These results are consistent with the spin-parity assignments of 0"'" for the _D3,j(2317) and 1+ 
for the L»sj(2457). 
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The narrow Ds'k'^ resonance at 2317 MeV/c^, recently 
observed by the BaBar collaboration Q , is naturally in- 
terpreted as a P-wave excitation of the cs system. The 
observation of a nearby and narrow D*Tr^ resonance by 
the CLEO collaboration Q supports this view, since 
the mass difference of the two observed states is con- 
sistent with the expected hyperfine splitting for a P- 
wave doublet with total light-quark angular momentum 
j = 1/2 0,0. The observed masses are, however, con- 
siderably lower than potential model predictions iQ] and 
similar to those of the cu j = 1/2 doublet states recently 
reported by Belle . This has led to speculation that the 
new Di*^Tr^ resonances, which we denote D^j, may be ex- 
otic mesons H H E [ll 111 El ■ Measurements of the 
Dgj quantum numbers and branching fractions (partic- 
ularly those for radiative decays) , will play an important 
role in determining the nature of these states. 

In this paper we report measurements of the D^j 
masses, widths and branching fractions using a sample of 
e+e" cc events collected with the Belle detector ^3 
at the KEKB coUider 

We reconstruct Df mesons using the decay chain 
Df — > (l)TT~^ and (j) K^K~ . To identify kaons or pions, 
we form a likelihood for each track, Ck{-k)i from dE/dx 
measurements in a 50-layer central drift chamber, the re- 
sponses from aerogel threshold Cerenkov counters, and 
time-of- flight scintillation counters. The kaon likelihood 
ratio, P{K/'k) — Ck/ICk + ^-r), has values between 
(likely to be a pion) and 1 (likely to be a kaon). 

For (j) K^K~ candidates we use oppositely-charged 
track pairs where one track has P^K/tt) > 0.5 and the 



other has P{K/tt) > 0.2, and with a K^K^ invariant 
mass that is within lOMeV/c^ (^ 2.50") of the nominal (f> 
mass. We define the (f> helicity angle 6h to be the angle 
between the direction of the K'^ and the D'^ in the cf) rest 
frame. For signal events this has a cos^ 9h distribution, 
while for background it is flat; we require | cos 9h\ > 0.35. 

We reconstruct Df candidates by combining a (f> 
candidate with a tt"*", which is a charged track with 
P{K/tt) < 0.9, and requiring Af((/)7r+) to be within 
lOMeV/c^ 2a) of the nominal _D+ mass. We use the 
D+ sideband regions 1920 < M(07r+) < 1940MeV/c2 
and 1998 < Af(07r+) < 2018MeV/c2 for background 
studies. 

For TT*^ reconstruction, we use photons with the e^e~ 
rest frame (CM) energies greater than 100 MeV and se- 
lect 77 pairs that have an invariant mass M(77) within 
10MeV/c^(^ 2(t) of the mass. For background 
studies we use the 7r° sideband regions 105 < M^^ < 
115MeV/c2 and 155 < M-y^ < 165MeV/c2. 

We reconstruct D*^ in the Df'j final state. We use 
photons with CM energies greater than 100 MeV and re- 
quire D*+ candidates to satisfy 127 < AM{D+j) < 
157MeV/c2(- 3(t), where AM{D+j) = M{D+-f) - 
MjjJr. The Z?*+ sideband regions are defined as 87 < 
AM{Dfj) < 117MeV/c2 and 167 < AM{Dfj) < 
197MeV/c2. 

The AM{D+TT°) = M{D+Tr°) - mass-difference 
distribution for Dfn'-* combinations with p*(Z3+7r°) > 
3.5GeV/c is shown in Fig. QJa). Here, and in analy- 
ses of other Dgj states and modes, we require the CM 
momentum to satisfy p*{Dsj) > 3.5GeV/c to remove 
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contributions from BB events. We do not remove mul- 
tiple candidates in the subsequent analysis. Also shown 
are the distributions for the Df (solid) and tt" (dashed) 
sideband regions. The prominent peak in the figure cor- 
responds to the Dsj(2317) Dfir'^ signal; the peak at 
small AM values is due to i:'*+(2112) ^ D+vrO. No peak 
is seen in the sideband distributions. 

Figure mb) shows the AAf(L>;+7rO) = M{D;+tt°) - 
M^,+ distribution for p*{D*+'K^) > 3.5GeV/c, where a 
peak corresponding to Dsj{2A57) i)*+7r° is evident. 
Also shown is the distribution for the -D*+ sideband re- 
gion, where we notice the presence of a wider peak in the 
Dsji2457) region. The AM{D;+Tr°) distributions for the 
Df and sideband regions show no such peak. 

To study the expected signal shape and detection ef- 
ficiencies, and determine the level of cross-feed between 
the two states, we use a Monte Carlo simulation that 
treats the Z?sj (2317) as a scalar particle with mass 
2317MeV/c2 decaying to Dfir^ and the i:)^j(2457) as 
an axial- vector particle with mass 2457MeV/c^ decay- 
ing to _D|+7r°. Zero intrinsic width is assigned to both 
states. We find that the Dsj{23l7) produces a peak of 
width 7.1 ± 0.2MeV/c^ in the AM{DfTr") distribution 
at its nominal mass, and a broader reflection peak (of 
width 12.3±1.8 MeV/c^) at a mass of 8 MeV/c^ above the 
Dsj{2457) peak. This latter peak corresponds to a Df 
and tt" from a Dsj{2317) decay that are combined with 
a random photon that passes the \M{D^'-f) — A/^*+| < 
15MeV/c^ requirement. (We refer to this as "feed- up 
background.") The I?sj (2457) produces a peak of width 
6.0 ±0.2 MeV/c^ at its nominal mass and a broader peak 
(of width 19.5 ± 3.6 MeV/c^), also at its nominal mass. 
The latter peak is due to events in which the photon 
from D*^ -D^7 is missed, and a random photon is 
reconstructed in its place (referred to as the "broken- 
signal background"). In addition, the I?sj(2457) pro- 
duces a reflection in the Dfn^ mass distribution with 
width 14.9±0.8MeV/c2 at a mass of 4MeV/c2 below the 
Dsj{2317) peak (referred to as "feed-down background" ) . 

While we must depend on the MC for separating 
the signal peak and the feed-down background in the 
Dsj{2317) region, the feed-up and broken-signal back- 
grounds for the Dsj{2457) region occur when £)*+7r" 
combinations are formed from candidates in the _D*+ 
mass sideband. This is evident in Fig.QIb). 

Figure EJb) shows the sideband-subtracted 



AM{D 



distribution together with the results 



of a fit that uses a Gaussian to represent the 15s,/ (2457) 
signal and a second-order polynomial for the background. 
The fit gives a signal yield of 126 ± 25 events with a 
peak value of AM = 344.1 ± 1.3MeV/c^ (corresponding 
to M = 2456.5 ± 1.3MeV/c2). The width from the 
fit, a = 5.8 ± 1.3 MeV/c^, is consistent with MC 
expectations for a zero intrinsic width particle. 

Figure |2Ia) shows the fit result for the Dsj{2317). 
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FIG. 1: (a) The AM{D+-k°) distribution. Data from the D+ 
(solid) and 7r° (dashed) sideband regions are also shown, (b) 
The AAI {Ds^ Tv°) distribution. Data from the Dl sideband 
(histogram) region is also shown. 
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FIG. 2: (a) The AM(D+7r°) distribution. The narrow 
Gaussian peak is the fitted Dsj {2317) signal, whereas the 
wider Gaussian peak is the feed-down background, (b) The 
AA4{Ds'^Tv^^) distribution after bin-by-bin subtraction of the 
D*'^ sideband from the D*^ signal distribution. The curve is 
the fit result. 



Here both the signal and the feed-down background are 
represented as Gaussian shapes modeled from the MC. 
The mean and a of the feed-down component are fixed 
according to the MC and normalized by the measured 
Dsj{2A57) yield. A third-order polynomial is used to 
represent the non-feed-down background. The fit gives a 
yield of 761 ± 44 events and a peak AM value of 348.7 ± 
0.5MeV/c2 (corresponding to M = 2317.2±0.5 MeV/c^). 
Here again, the width from the fit, a = 7.6±0.5 MeV/c^, 
is consistent with MC expectations for a zero intrinsic 
width particle. 

There are systematic errors in the measurements due 
to uncertainties in the: i) tt'^ energy calibration; ii) pa- 
rameterization of the cross-feed backgrounds; iii) param- 
eterization for the non-cross-feed backgrounds; iv) possi- 
ble discrepancies between the input and output seen in 
the MC simulations; and v) the uncertainty in the world 
average value for M^+ and Mjj,+ . 

The tt'^ energy calibration is studied using 
£'*+(2112) ^ -Dj^Tr" events in the same data sam- 
ple. We measure AM = 144.3 ± O.lMeV/c^ and 
a = 1.0 ± O.lMeV/c^, which agrees well with the 
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PDG value of AM = 143.8 ± 0.4MeV/c2. 
MC, which uses the PDG value as an input, 



The 
gives 



AM = 143.9 ± 0.1MeV/c2 and ct = 1.0 ± O.lMeV/c^. 
(The errors quoted here are statistical only). We 
attribute the difference to the 7r° energy calibration 
uncertainty, and conservatively assign a ±0.6MeV/c^ 
error to this effect. This error only contributes to the 
mass measurements. 

For the cross-feed background to the Dsj{23l7) signal, 
we vary the feed-down background parameters and the 
Dsj(2457) yield by ±lcr and assign the variation in out- 
put values as errors. For the Z)sj(2457), we determine 
the uncertainty of the feed-up fraction from the differ- 
ence between the D* signal region and the sideband re- 
gion using the MC. For the non-cross-feed background, 
we repeat the fit using a second-order polynomial for the 
Dsj{2317) and a linear function for the I?sj(2457) and 
assign the difference as errors. Differences between the 
MC input and output values for the Dsj parameters can 
reflect possible errors arising from the choice of signal 
shape and other factors in the analysis. We assign these 
differences as errors. 

The final results for the masses are 



M(i^,j(2317)) 
M{Dsj{2A57)) 



2317.2 ± 0.5(stat) ± 0.9(syst) MeV/c^ 
2456.5 ± 1.3(stat) ± 1.3(syst) MeV/c^. 



The M(£)sj(2317)) result is consistent with BaBar Ij 
and CLEO results 0. Our M(i:>,j(2457)) value is 
consistent with BaBar ^3 but significantly lower than 
that from CLEO 0. We set upper limits for the 
natural widths of r(£)^j(2317)) < 4.6MeV/c2 and 
r(L>,j(2457)) < 5.5MeV/c2 (90% C.L.), respectively. 

Using the observed signal yields of 761 ± 44(stat) ± 
30(syst) and 126 ± 25(stat) ± 12(syst)for the L>sj(2317) 
and Dsj(2457), and the detection efficiencies of 8.2% and 
4.7% for the Ds.7(2317) and Dsj{2457), we determine the 
ratio 

a(i?,j (2457)) • S(i?+ (2457) ^ D*+7r°) 



a{D s J {2317)) ■B{D+j{2317)- 
0.29 ± 0.06(stat) ± 0.03(syst). 



D+nO) 



The detection efficiencies are determined from the MC 
assuming the same fragmentation function for the two 
states. The dominant source of systematic error is the 
systematic uncertainty in the Dsj{2A57) yield. 

In the Ds,7(2457) region of the AM(D+7r°) distribu- 
tion, we find 22 ± 22 events from a fit to a possible 
-Ds J (2457) signal. From this, we obtain the upper limit 



B(i^+ (2457) ^ D+TT^) 



BiD 



."7(2457) 



Dt+ir") 



< 0.21 (90% C.L.). 



The decay to a pseudo-scalar pair is allowed for a state 
with a parity of (—1)"^. Thus, absence of such a decay 
disfavors Ds,/(2457) having of 0+ or 1". 
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FIG. 3: (a) The AM(D+7) distribution. The curve is a 
fit using a double Gaussian for the signal and a third-order 
polynomial for the background, (b) The AM(_D^7r'''7r~) dis- 
tribution. The curve is a fit using Gaussian for the signals 
and a third-order polynomial for the background. 



Figure Efa) shows the AM{Df-f) M{Dfj) - M^+ 
distribution. Here photons are required to have ener- 
gies greater than 600 MeV in the CM and those that 
form a tt" when combined with another photon in the 
event are not used. A clear peak near AAI{Df'-f) ~ 
490MeV/c^, corresponding to the Dsj{2A57), is ob- 
served. No peak is found in the Dsj{2317) region. The 
Df sideband distribution, shown as a histogram, shows 
no structure. We fit the distribution with a double 
Gaussian for the signal, which is determined from the 
MC, and a third-order polynomial for the background. 
The fit yields 152 ± 18 (stat) events and a AM peak 
at 491.0 ± 1.3(stat) ± 1.9(syst) MeV/c^ (corresponding 
to M = 2459.5 ± 1.3(stat) ± 2.0(syst) MeV/c^). The 
-Dsj(2457) mass determined here is consistent with the 
value determined from D*tt'^ decays. 

Using the detection efficiency of 10.2% for the Dfj 
decay mode, we determine the branching fraction ratio 



g(i^+ (2457) ^ J+7) 
S(i?+ (2457) ^ DS+ttO) 



= 0.55±0.13(stat)±0.G8(syst). 



This result, which has a statistical significance of lOcr, 
is consistent with the first measurement by Belle [l^ 
0.38 ± O.ll(stat) ± 0.04(syst) with B DD^{2A57) de- 
cays, and with the theoretical predictions [3,[l3- The 
existence of the Dsj{24:57) Z3s7 mode rules out the 0^ 
quantum number assignments for the D^j (2A57) state. 
For the Dsj{2317), we obtain the upper limit 



g(i?+ (2317) ^ D+j) 
B{D+,{2in) ^ i?,7r") 



< 0.05 (90% C.L. 



From the M{D*+-i) = M{D*+"f) 
we determine the upper limits 



M 



distribution, 



B{Dtj{2il7) ^ D*+^) 
B{D+j{2il7) ^ D,TT^) 



< 0.18 (90% C.L.) and 



B{Dtj{2457) 



S(i?+ (2457) ^ Dt 



< 0.31 (90% C.L.). 
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FigureEIb) shows the AM(£)+7r+7r-) = M(D+7r+7r-) - 
Mjj+ distribution. For pions, we require at least one 
of them to have momentum greater than 300MeV/c 
in the CM, one with P{K/tt) < 0.1 and other with 



P{K/Tr) < 0.9, and |M(7r^ 



Mks\ > 15MeV/c2 



A clear peak near AM(D+7r+7r") - 490MeV/c2, corre- 
sponding to the Dsj{24:57), is observed. Evidence of an 
additional peak near AM(D+7r+7r~) ~ 570MeV/c^ cor- 
responding to Dsi(2536) is also visible. No peak is found 
in the Dsj{2317) region. The sideband distribution, 
shown as a histogram, shows no structure. We fit the dis- 
tribution with Gaussians for the signals, which are deter- 
mined from the MC, and a third-order polynomial for the 
background. The fit yields 59.7 ± 11.5(stat) events and a 
AM peak at 491.4 ± 0.9(stat) ± 1.5(syst) MeV/c^ (corre- 
sponding to M = 2459.9 ± 0.9(stat) ± 1.6(syst) MeV/c^) 
for Dsj{2A57), and 56.5±13.4(stat) events for Dsi{2536). 
The statistical significance is 5.7cr for Dsj{2A57), and 
4.5(7 for Z?si(2536). This is the first observation of the 
Dsj{2i57) ^ D+TT+TT" decay mode. 

The existence of the I?sj(2457) DsTt^tt^ mode also 
rules out the 0+ assignment for Dsj{24:57). Using the de- 
tection efficiency of 15.8% for the D^Tr+Tr" decay mode, 
we determine the branching fraction ratio 



g(ZJ+(2457) 



(2457) -> Dt+nO) 



0.14±0.04(stat)±0.02(syst), 



where the systematic error is dominated by the system- 
atic uncertainty of the Z3sj(2457) Z3*+7r'^ yield. We 
establish the upper limit 



6(i?+ (2317) ^ Di 



BiD+j{2317) ^ i^+7rO) 



^ < 4 X 10"^ 



(90% C.L.). 



Using the detection efficiency of 14.3% for the 
Dsi{2536) Ds7r+7r~ decay mode which assumes the 
same fragmentation function for the _Dsi(2536) and 
-Dsj(2457), we establish the cross section times branching 
fraction ratio 

cr(L>^i(2536)) • 6(15+1 (2536) ^ D+tt+t:-) 



a{Dsj{2i57))-BiD+j{2457)- 
1.05 ± 0.32(stat) ± 0.06(syst). 



In summary, we observe radiative and dipion decays of 
the -Dsj(2457) and set upper limits on the corresponding 
decays of the i:»^j(2317). We determine the i:>s,/(2317) 
and Dsj{2A57) masses from their decays to D^tt'^ and 
Z)*+7r°, respectively, and set an upper limit on the de- 
cay of _Dsj(2457) to -D+tt". These results are consistent 



with the spin-parity assignments for the I?s,7(2317) and 
Dsj{2457) of 0+ and 1+, respectively. 
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